Source Code for In-person schooling and associated COVID-19 risk in the United States over Spring Semester 2021
Kirsten E. Wiens, Claire P. Smith, Elena Badillo-Goicoechea, Kyra H. Grantz, M. Kate Grabowski, Andrew S. Azman, Elizabeth A. Stuart, Justin Lessler
Code is provided to reproduce the primary in-person schooling (1-Main-Schooling-Analysis.Rmd) and occupation/work outside home (2-Educators-Analysis.Rmd) analyses using a synthetic dataset. Hence, results produced when knitting these documents will not match those presented in the accompanying paper.
Data to reproduce paper figures/tables are freely available from the CMU Delphi Research Group to researchers at universities and non-profits as detailed at Getting Data Access - Delphi Epidata API (cmu-delphi.github.io).
Please reach out to justin[at]jhu.edu if you have any questions.
Table S1
Sociodemographic characteristics of participants with ≥1 school-aged child in the household, comparing those reporting no in-person schooling to those reporting any part-time and full-time in-person schooling. Observed and survey-weighted percentages reported.
Table S2
Monthly COVID-19-related outcomes in different subsets of US CTIS respondents, including total number, observed percent, and percent weighted to account for survey design.
Figure 1
Changes over time in in-person schooling by county. Distribution of survey responses from Jan. 12 to Mar. 31 (left column), Apr. 1 to Jun. 12 (center column) and change over time (right column). Results are shown for (A) number of survey respondents reporting ≧1 school-aged child in the household, (B) percent reporting in-person schooling, (C) percent of respondents with in-person schooling reporting full-time in-person instruction, and (D) average number of school-based mitigation measures. “Relative Amount” in the right column indicates values in Apr. 1 to Jun. 12 (center column) divided by values in Jan. 12 to Mar. 31 (left column). Grey indicates county/periods where fewer than 10 respondents reported in-person schooling.
Figure 2
Concurrent changes in in-person schooling, vaccination, and variant prevalence. Changes over time between Jan. 12 and Jun. 12 in (A) percent of respondents living with school-aged children reporting any in-person schooling, (B) average number of school-based mitigation measures, (C) smoothed percent of GISAID SARS-CoV-2 sequenced isolates that were Alpha or Delta, and (D) percent of respondents living with school-aged children reporting having received at least one COVID-19 vaccine dose. Upper panel line plots show national averages by week weighted using US CTIS survey weights (A-B), state population (C), and county population (D). Lower panel histograms show the number of county-months with the indicated percentages (A,C,D) or numbers (B). Note that the number of respondents decreased over time (Fig 1d), which may contribute to the increasing number of zero values in the histograms in June.
Movie S1
Changes in in-person schooling, vaccination, and variant prevalence. Top left: County-level smoothed weekly estimates of percent of individuals living with school-aged children participating in part-time or full-time in-person schooling. Top right: County-level smoothed weekly estimates of the percent of individuals living with school-aged children that had received any number of COVID-19 vaccine doses. Bottom left: State-level smoothed weekly estimates of the proportion of circulating SARS-CoV-2 strains that were the Alpha variant. Bottom right: State-level smoothed weekly estimates of the proportion of circulating SARS-CoV-2 strains that were the Delta variant. In-person schooling and vaccination estimated using United States COVID-19 Trends and Impact Survey data; counties with fewer than 20 respondents reporting in-person schooling excluded. Variant prevalence estimated using GISAID sequencing data. Dates represent the last day of each week.
